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Abstract 

It is shown that the physical constraint of the Anomalous Green's function 
gives a natural pairing condition. The resulting self-consistency equation is 
directly related to the BCS gap equation. Both inhomogeneous and homoge- 
neous systems are considered to illustrate the importance of the constraint. 
Especially we find weak localization correction to the phonon-mediated inter- 
action. 
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1. Introduction 



Recently Green's function treatment of impurity effects on the superconductors was 
shown to differ from that of the BCS theory. ^'^'^ The electron density of states change caused 
by nonmagnetic impurities is overestimated in the Abrikosov and Gor'kov's (AG) theory,^ 
which is inconsistent with Anderson's theorem.^ (Details of the strong coupling theory of 
impure superconductors are discussed elsewhere.^) For magnetic impurity effects, Kim and 
Overhauser^ proposed a BCS type theory with different predictions: (i) The initial slope of 
Tc decrease depends on the superconductor and is not the universal constant.^ (ii) The Tc 
reduction by exchange scattering is partially suppressed by potential scattering when the 
overall mean free path is smaller than the coherence length. This compensation has been 
confirmed in several experiments. '^■^'^■^^'^^ 

In this letter we show that the physical constraint of the Anomalous Green's function 
gives a natural pairing condition in Gor'kov's formalism. The resulting self-consistency 
equation is nothing but another form of the BCS gap equation. Accordingly the above 
discrepancies are settled in favor of the BCS theory. 

To illustrate, we first consider inhomogeneous systems with nonmagnetic and magnetic 
impurities. Weak localization correction to the phonon-mediated interaction is derived. We 
also discuss a homogeneous system. 

2. Inhomogeneous System: Nonmagnetic Impurity Case 

Let's consider an electron gas in the presence of the ordinary impurities. In the second- 
quantized representation, the Hamiltonian is 




where Uo{y) — Uo5{v — Rj) is the impurity potential. The field operator, '^{va)^ may be 
expanded by the exact scattered states basis set '?/'n(r), that is, 

^(ra) = ^^/;n(r)cna, (2) 

n 
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where Cna is a destruction operator of the electron. The Gor'kov's Anomalous Green's 
function F(r, r') is defined^^ 

F(r,r') = -<^(rT)^(r'i)>- (3) 
Accordingly, the pair potential is given 

A(r) = 1/F(r,r) = V"^F(r,r,u;). (4) 

The a;'s are ujn — (2n + l)7rT for all integer n. Near the transition temperature, it has been 
believed that^"^ 

F(r, r', uj) = l A(l)Gi(r, l)Gi,(r', l)dl, (5) 



where 



Gt(r,l) = sM£>i;i(l), (6) 



and 



luj — e„ 



Gi„(r',l) = i:*^^^^. (7) 

Substituting Eq. (5) in Eq. (4), we find a well-known self- consistency equation 

A(r) =VTY.f A(l)GL(r, l)Gi^(r, \)d\. (8) 

2.1 A Constraint on the Anomalous Green's function 

Now we consider the physical constraint of the Anomalous Green's function. If we 
average over the impurity positions, the system becomes homogeneous. Consequently, the 
Anomalous Green's function should be a function of r — r' after the impurity average, i.e.. 



, r^mp —. -primp , . 

F(r, r'j = F{r - r') . (9) 
—imp jngans an average over impurity positions Rj. 



However, if we substitute Eqs. (6) and (7) into Eq. (5) we find extra pairings between 
n t and n'{j^ n) |, which violate the physical constraint of the anomalous Green's function. 
n denotes the time reversed partner of the scattered state n. Using the scattered states 

V'„w(r) = e^^-' + E -^^jE e-^^-^'le'M-- + • • • , (10) 

q ^k+q i 

it is easy to show that 

q,i y^k H+q)y^k' '^k'-q) 

^f{r-r'), (11) 

and 



k ~ ^k+q' 

= g(r-T'). (12) 



q,i \^k k+q) 



Note that even the zeroth order term of extra pairings between n \ and n'{^ n) [ is not a 
function of r — r'. Prom Eqs. (5) and (11) one finds that 

F{r, r', u) ~ V'„(fe)T(r)V'„'(ifcOi(r') 



^F{r-r',ojr', (13) 

in the presence of the extra pairing terms. 

It is clear that only Anderson's pairing between n | and n J, is compatible with the 
homogeneity condition of the Anomalous Green's function after averaging out the impurity 
positions. Accordingly, the physical constraint of the Anomalous Green's function gives 
a natural pairing condition. By eliminating the extra pairing, the 'corrected' Anomalous 
Green's function is given 

F(r, r', 00)^1 A(l){Gl(r, l)GiUr', l)},.,d\, (14) 

where 



{G'jr, l)GL„(r', 1)}„ = E E *^i^^^^*„..„ (16) 

and p.p. means proper pairing constraint, which dictates pairing between n t and n J,. 
Consequently, the self-consistency equation also needs a proper pairing constraint derived 
from the Anomalous Green's function. The resulting equation is 

A(r) ^VTY^f A(l){Gi(r, l)Gi,(r, 1)},.,^!. (16) 

Let's compare Eq. (16) with the BCS gap equation, 

K--Y.ynn'^tanhC^l (17) 

n' 

where^'^^ 

Kn' ^-V j rn'ir)rMM^)M^)dr. (18) 
In fact, the relation between A(r) and A„ was obtained by Ma and Lee,^^ (i-e.), 

A„ = I <(r)V^*(r)A(r)dr. (19) 
Substitution of Eq. (16) into Eq. (19) leads to the BCS gap equation Eq. (17), since^^ 

Note that because Eqs. (16) and (17) are basically equivalent, they give the same result 
for the impurity effect on the superconductors. Up to the first order in the impurity con- 
centration, one finds^'^ Vnn' — —V and the density of states change is 0{1/Eft). Ep is the 
Fermi energy and r is the scattering time. Therefore Anderson's theorem is obtained. In 
the AG theory^ the density of states change is 0{1/ujdt),^ which gives a large decrease of 
Tc- is the Debye energy. The Dyson expansion of Green's function in the presence of im- 
purities can not distinguish two energy scales ujd and Bp coming from the phonon-mediated 
interaction and the virtual impurity scattering. 

2.2 Weak localization correction 



We are now ready to consider weak localization^''''^^'^^ correction to the phonon-mediated 
interaction. Mott and Kaveh^^ showed that the wavefunctions for the weakly locahzed states 
may be written as a mixture of power-law and extended wavefunctions, of the form, for three 
dimensions, 

V'f(r)=>l3*L + 53^, (21) 



and, for two dimensions. 



where 



^^''(r) = A2¥„, + B2^, (22) 



^l-l-4.i3|(i4). B,-i^,, (23) 

and 

Al = l-27TBlln{L/i), B2^^r-o- (2^) 

TT K jpZ 

£ and L are the elastic and inelastic mean free paths. 

We can find the effective interaction by substituting Eqs. (21) and (22) into Eq. (18). 
However, we must be careful in calculating contributions from the power-law wavefunction. 
Because we are concerned with the bound state of Cooper pairs in a BCS condensate, only 
the power-law wavefunctions within the the BCS coherence length are relevant.^ This 
is analogous to the insensitivity of the localized states with the change of the boundary 
conditions. Accordingly, it is obtained 

V^. = -VlAl + (|)8x^|i3|(l _ i) + 0(Bieja?)] 



Note that ^^/^ Co/^ second terms of the right hand side of Eqs. (25) and (26) 

denote the probabihties of finding the centers of the power-law wavefunctions within the 
Cooper pair radius in 3 and 2 dimensions, fl is the volume of the system. It is remarkable 
that the same correction occurs in the conductivity and the phonon-mediated interaction 
by weak localization. Consequently, one expects in 1 dimension 

where a is the radius of the wire. 

There are many experimental results which show the reduction of Tc caused by weak 
localization. ^^'^^'^'^ Previously, it was interpreted by the enhanced Coulomb repulsion. How- 
ever, Dynes et al.^^ found a decrease of the Coulomb pseudo-potential /i* with decreasing Tc. 
We believe that this signals the importance of weak localization correction to the phonon- 
mediated interaction. Comparisons with the experiments will be published elsewhere. 

3. Inhomogeneous System: Magnetic Impurity Case 

Wc now consider the effect of magnetic impurities on the superconductivity. Magnetic 
interaction is given by 

Um{r) = ^Js-S,S{r-Ri), (28) 

i 

where s — \a. The three components of a are the Pauli matrices. For convenience, we 
consider only the z-component of t/rn(r). The effect of each x and y-component is basically 
the same with that of z-component.^ Notice that up and down spin electrons feel different 
potential because of the Pauli spin matrix, cr^, i.e., 

+ E ^i^Scose5{v - mrn^{v) = KtCtlr), (29) 
- E pcose5{v - R.)]<i(r) = S„i<i(r), (30) 



where S = yS{S + 1) and 6 is the polar angle of the fixed local-spin S. As a result, one 
finds, 



The Anomalous Green's function, near the transition temperature, is given by^ 
F(r,r',u;) = / A(l)G'I(r, l)G'i,(r', 

where 
and 



As in Sec. 2.1, we check the physical constraint of the Anomalous Green's function 



- rimp r^mp 

F{y,r',u) ^ F[y — y',uj) 
Using the scattered states 

and 

q ^k+q i 

it is straight forward to show that 



g«<?-(r-r') 



iJ%~i^W^~J^"''' = e^('="+'='-'-'Hl - {JScose/2f E 7 ^T7 7 + 

7^/(r-r'), 



and 



The minus sign in the perturbation correction term of Eqs. (37) and (39) is the noteworthy 
feature of exchange scattering. Prom Eqs. (35), (38) and (39), it is obvious that only the 
pairing between the scattered states n{k) j and n{—k) I is allowed. Note that n{k) t and 

— * 

n{—k) I are degenerate, 

E-^-" = Err". (40) 

The resulting self-consistency equation is 

A(r) ^VTY^j A(l){Gi(r, l)Gi,(r, 1)^^!. (41) 



where 



{Gl(r, l)Gi.(r, 1) U. = E f^Wil^ E ''''^^'-"^^ 6r,.,. (42) 



Accordingly, we obtain the BCS gap equation 



A„ = -EKn'^tonM^), (43) 



where 



By taking an impurity average of Eq. (44), we find that 

Vnn' = -V[l-2Zi^ _ -2Zi}^——- -J, (45) 

q y'k ^k+q> q V'k' ^k'+q) 

where is the number of the magnetic impurities. Note that this result is the same as the 
z-component in Eq. (55) of Kim and Overhauser.^ In other words, eliminating the extra 
pairing terms from the AG theory^ gives rise to Kim and Overhauser's result.^ 

4. Homogeneous System 

Near the transition temperature, the Anomalous Green's function is given by 

F(r,r',a;) = j A(l)GL(r, l)Gi,(r', l)dl, (46) 



where 



0r,(r)0i (1) 



and 



k' 



-luj - 



Note that 0g(r) = e^^'^ . Accordingly, the self-consistency equation is 

A(r) ^VTJ2J A(l)GT(r,l)Gi,(r,l)dl 



K{r, 1) A(l)dl. (49) 

For a homogeneous system, Gor'kov^^ pointed out correctly that F{r, r') should depend 
only on r — r', that is, 

F(r,r') = F(r-r'). (50) 

However, observe that 

</'.-T(r)</>.-u(r') = e^^^^-^+^"'-'\ (51) 

and 

Eqs. (46) and (49) include the extra pairing terms between k t and k' | (7^ — /c j), which 
do not satisfy the homogeneity condition of Eq. (50). In this case, the self-consistency 
condition of the pair potential happens to eliminate the extra pairing in Eq. (49) because 
of the orthogonality of the wavef unctions. 

However, it is important to eliminate the extra pairing in the Anomalous Green's function 

— * 

from the beginning. Note that the kernel K{r, 1) is not for the pairing between k t and 
—k I, but for the pairing between the states which are the linear combination of the plane 
wave states 0g(r).^^ From the Bogoliubov-de Gennes equations, it can be shown that the 
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corresponding unitary transformation leads to the vacuum state where the pairing occurs 
between the states, which are the hnear combination of the plane wave states. The proper 
kernel for the BCS pairing between k j and —k I is 

K'^ir, 1) = VT E{Gl(r, IjGi^r', l)}p.p., (53) 

where 

k k' 

Consequently, the 'corrected' Anomalous Green's function is 

F(r,r',^) = j A(l){GL(r,l)Gi,(r',l)}^.^.cil. (55) 

In fact, Eq. (55) is not new. Another form of this equation was already obtained as 
shown^''"'^^ 

F(r, r') = ^tanh{^)exp[ik ■ (r - r')]. (56) 



2er ^2T' 
k ^ 



By using 



^k = / 0;-(r)r_g(r)A(r)dr, (57) 

and 

(58) 

Eq. (56) is transformed to 

F(r,rO=Enr,r',u;), 

= J2[ A(l){Gl(r,l)Gi,(r',l)U.a (59) 

Therefore, Eq. (55) is confirmed. 

5. Conclusion 
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It is shown that the physical constraint of the Anomalous Green's function gives a natural 
pairing constraint. The resulting self-consistency equation is nothing but another form of the 
BCS gap equation. Anderson's pairing between the time reversed scattered state partners 
are obtained in the presence of the ordinary impurities. Weak localization correction to the 
phonon-mediated interaction is calculated. In the presence of the magnetic impurities, the 
degenerate partners are paired, which give rise to the result of Kim and Overhauser. In the 
case of a homogeneous system, the BCS pairing is required. 
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